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I. INTRODUCTION 

Experimental and theoretical studies of spin polarized 
transport in hybrid magnetic nanostructures is a fron- 
tier in mesoscopic physics. The most prominent example 
of conceptual, technological, and commercial impact is 
the giant magnetoresistance effect utilized in magnetic 
information storage devices. In order to gain a deeper 
understanding of spin and charge transport, and to en- 
hance circuit functionality and efficiency, more complex 
structures are fabricated and studied. In recent years, 
hybrid nanoscale circuits containing normal conductors, 
ferromagnets, and superconductors have been realized. 
These structures allow observation and understanding of 
competing mechanisms of electron-electron interactions. 

The simultaneous existence of ferromagnetism and su- 
perconductivity is rare. In ferromagnets, the exchange 
interaction lifts the spin-degeneracy and induces an itin- 
erant spin polarization. In s-wave superconductors, on 
the other hand, electrons with anti-parallel spins form 
Cooper pairs. In homogenous conventional ferromag- 
nets (Fe, Ni, Co, and alloys thereof), the large exchange 
interaction efficiently dephases electron-hole pairs, and 
eliminates singlet superconducting correlations over dis- 
tances larger than the ferromagnetic coherence length. 
This would suggest a short-range superconducting prox- 
imity effect in transition metal ferromagnets Such a 
simple picture cannot explain recent measurements on Co 
and Ni ferromagnets coupled to Al superconductors, how- 
ever, where a substantial resistance drop was observed at 
the onset of superconductivity^ The simple picture also 
fails to explain the long-range superconducting proximity 
effect recently observed via the Josephson supercurrent 



through a half-metallic ferromagnet Subsequent the- 
oretical work show that induced triplet superconducting 
correlations give rise to long ranged proximity effect in 
transition metal ferromagnets £2 Triplet superconduct- 
ing correlations are insensitive to the pair-breaking ex- 
change interaction and exhibit a longer coherence length, 
similar to that of superconducting correlations in nor- 
mal metals. It is now established that spin-flip pro- 
cesses in a ferromagnet can convert singlet into triplet 
pair correlations. A spatially inhomogeneous magnetiza- 
tion texture^ or magnonsii~— are examples of spin-flip 
sources that are able to induce long ranged triplet corre- 
lations. 

In this report, we focus our attention on the influ- 
ence of magnons on the transport properties in normal 
metal-ferromagnet-superconductor systems. Even nor- 
mal metal-ferromagnet systems without superconductors 
exhibit intriguing physics, and especially the interaction 
between spin and charge currents and the magnetic or- 
der parameter in such structures have attracted tremen- 
dous interest. For instance, a non-collinear spin flow to- 
wards a ferromagnet exerts a torque on the magnetiza- 
tion, a spin transfer torque, that can excite the magneti- 
zation and even induce steady state, precessional motion 
of the ferromagnetic order parameterj 14 ' 15 The inverse 
effect is also of significant interest: A precessing ferro- 
magnet in electrochemical equilibrium with its environ- 
ment, acts as a "spin battery" by emitting (or "pumping") 
pure spin currents into neighboring materials.-^ When 
emitted spins are dissipated in adjacent materials, spin 
pumping enhances magnetic dissipation in the precess- 
ing ferromagnet, and thus increases observed linewidths 
in FMR experiments^ 

Some ideas from spin transfer physics in normal 
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metal-ferromagnet structures were recently used to 
study superconductor-ferromagnet systems. A FMR 
experiment and the following theoretical analysis 19 
have shown how spin pumping can be used to visualize 
proximity effects and spin relaxation processes inside the 
superconductor. In essence, in metallic contacts, ferro- 
magnetic correlations reduce the superconducting order 
parameter close to the layer interface, enabling pumped 
sub-gap electrons to enter and deposit spin in the super- 
conductor. This is a prime example of how the inverse 
proximity effect affects the FMR linewidth broadening 
when typical spin-flip lengths are comparable to the su- 
perconducting coherence length. 20 

We direct our attention to a different aspect of the in- 
terplay between magnetization and carrier dynamics in 
ferromagnet-superconductor structures. In contrast to 
the works mentioned above, where the magnetization dy- 
namics have been the primary concern, we will consider 
how a precessing magnetization and an applied voltage 
bias induce spin and charge currents in a normal metal- 
ferromagnet-superconductor (N|F|S) trilayer. The com- 
puted charge currents can be measured directly, whereas 
spin currents can possibly be measured by its dissipa- 
tive effect on the precessing ferromagnet, its spin trans- 
fer torque effect on a second ferromagnet, or via spin- 
filtering as a charge buildup on another ferromagnet.— 
Related to our work, sub-gap transport properties have 
recently been studied in a normal metal-ferromagnetic 
superconductor structure. 21 In ferromagnetic supercon- 
ductors, magnetic and electron-hole correlations coex- 
ist which can result in novel transport and dynami- 
cal magnetic phenomena. It was shown how supercon- 
ducting correlations, namely Andreev reflections at the 
layer interface, add features to the results of spin and 
charge pumping in normal metal-ferromagnet systems. 
In this report, we also consider how pumping in the 
N|F|S trilayer is related to pumping in the normal metal- 
ferromagnetic superconductor system as studied in Ref. 

m 

Diffusive transport in hybrid superconductor-normal 
metal systems is usually formulated within a quasiclassi- 
cal description* 2 ^ Although this description give qualita- 
tive insight into transport properties of superconductor- 
ferromagnet systems ) 10 ! 13 the formalism is limited to fer- 
romagnets with exchange interactions much smaller than 
the Fermi energy. Thus, a quasiclassical description can- 
not be used to quantitatively study transport in transi- 
tion metal ferromagnets Fe, Ni and Co used in exper- 
iments. This is one of the reaons why we adopt the 
scattering theory to transport^ 3 - Another reasone is that 
scattering theory captures adiabatic slow time-dependent 
variations of the magnetization direction well. 

Scattering theory has proven most useful in the study 
of stationary charge and spin currents in magnetoelec- 
tronic structures^! and the time-dependent generaliza- 
tion has successfully been applied to describe para- 
metric pumping of charge 2 ^ - — and spin currents^ 
For the N|F|S structure under consideration, we de- 



rive charge- and spin currents in the normal metal 
conductor in response to a slowly precessing ferromag- 
netic exchange field and applied bias voltage. We focus 
on sub-gap energies, and how Andreev scattering con- 
tributes to the conductivites of the currents. In electro- 
chemical equilibrium, we make contact with the results 
for pumping in normal metal-ferromagnetic superconduc- 
tor structures? 2 ! We proceed by detailing how time- and 
energy gradients of the total scattering matrix contribute 
to non-equilibrium pumped currents, and find that both 
charge and longitudinal spin currents are unaffected by 
the precessing magnetization. Finally, we consider non- 
equilibrium charge and spin currents for trilayers where 
the ferromagnetic region is longer than the transverse 
spin coherence length. 

This paper is organized in the following way: The 
N|F|S system is described in Sec. HI1 In Sec. we 
use time-dependent scattering theory to derive general 
expressions for charge and spin currents to first order 
in pumping frequency. The total scattering matrix for 
the system is then invoked in Sec. IIVI to obtain non- 
equilibrium pumped currents. Our conclusions are in 
Sec. H 



II. MODEL DESCRIPTION 

The system is sketched in Fig. [1] It consists of a su- 
perconductor (S) in series with a ferromagnet (F) and 
a normal metal lead (Ni). Ni is ideally coupled to a 
normal metal reservoir (N res ). We assume N ros and S 
to be in local thermal equilibrium, and denote a pos- 
sible chemical potential difference between the normal 
and the superconducting side as /ijv — /is = eV. Spin- 
orbit interactions are disregarded, and the ferromag- 
netic order parameter is assumed to be homogeneous 
and with a fixed magnitude A xc inside F. Its direc- 
tion is along the time-dependent unit vector m(t) — 
(sin#(i) cosfii, sin9(t) sinQi, cos0(t)). The precessing 
magnetization serves as the pumping parameter in the 
system. 
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FIG. 1: A ferromagnetic scattering region (F) is connected to 
a superconductor (S) and a normal metal reservoir (N rcs ) via 
two normal metal leads (Ni and Na). Amplitudes of outgoing 
(incoming) carrier states are given by b (a). 
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We focus on sub-gap transport properties. Thus, pos- 
sible scattering processes include Andreev reflections at 
the F|S interface 2 ^, and spin-dependent normal scattering 
inside F. Following a standard procedure^ the scattering 
problem is greatly simplified by utilizing spatially sepa- 
rated regions where scattering processes occur. This is 
achieved by inserting a fictitious normal metal lead (N2) 
between F and S. We assume that N2 is longer than the 
Fermi wavelength, so that asymptotic, plane wave solu- 
tions are applicable in this region. The total scattering 
matrix is a concatenation of the scattering matrices for 
Ni |F|N2 and for Andreev reflections at the N2 |S interface. 
Transport between F and S is mediated by the ballistic 
N 2 lead. 

The singlet superconductor is described by the BCS 
Hamiltonian 

U= / dr*t(r)ff (r)* CT (r) 

<r=t,,r 

+ fdr {A(r)*}(r)*|(r) + A*(r)*4,(r)*t(r)}, (1) 

where Ho is the normal state, single-particle Hamiltonian 
and A(r) the superconducting gap. We model the gap 
by a step function, A(r) = Ae 1 *9(x), where the phase 
(j> is constant, and x is the coordinate perpendicular to 
the N2IS interface. We take the Fermi energy Ep to 
be the largest energy scale, and focus on the low energy 
transport properties in regimes when the superconduct- 
ing gap is much less than the exchange interaction in 
the ferromagnet A xc , eV < A <C A xc ,Ep. The Hamil- 
tonian (QJ is diagonalized by the following Bogoliubov 
transformation^ 

iW =^{%Un(r,*)+nr*v*(r,a)} , (2) 



where 7„ are quasiparticle annihilation (creation) oper- 
ators that satisfy the fermionic anti-commutation rela- 
tion 



III. TIME-DEPENDENT SCATTERING 
THEORY 

We now focus on the time-dependent scattering the- 
ory for the N|F|S structure in Fig. JH apply the general 
framework established in Refs. [ll2l-|2V|l3J, and make 
use of the scattering theory for hybrid superconductor- 
normal metal structures discussed in Rcfs. 29.32, We find 
it most convenient to study a slowly precessing magne- 
tization by a scattering matrix expressed in the Wigner 
representation, 22 making the derivation of pumped cur- 
rents similar to that carried out for normal systems in 

Refs. mm. 

In order to describe a scattering potential of arbitrary 
time-dependence, we start by considering the two-time 
scattering matrix S(t,t'), that relates annihilation op- 
erators between states outgoing and incoming from the 
scattering region: 



(5) 



As indicated in Fig. [TJ b a : (a a ) annihilates the outgoing 
(incoming) state a. a labels electron-hole Nambu space 
index, spin and transverse wave-guide number. We as- 
sume that the reservoirs connected to the scattering re- 
gion are in local thermal equilibrium, and that incoming 
carriers from the normal metal reservoir fulfill 



(al(e)a a '(e')) = 8 aya ,5(e - e')f a (e), 



(6) 



where the brackets indicate a quantum and statistical 
average, and 



f e(h) (e) = f (e-a e ^eV) 



1 + e 



(e-cr° {h) eV)/k B T cl 



> (7) 



{7m, 7^} = Sr, 



(3) 



where cr e W = +(— )1, and f e (h)( £ ) is the Fermi-Dirac dis- 
tribution of incoming electrons (holes) at a charge bias 
eV and electron temperature T e ; . We will eventually con- 
sider electron temperature to be lower than the super- 
conducting gap. the We will now proceed by computing 
charge and spin currents in the system. 



The transformation @ results in a matrix equation for 
the quasiparticle eigenfunctions u n and v n : 



H (r) iA(r)a y 
-iA* (r)a y -H$(r) 



Un(r) 

v n (r) 



u n (r) 
v n (r) 



(4) 



The quasiparticle excitation energy s n is measured with 
respect to the chemical potential of the superconductor, 
which is set to zero. a y is a Pauli matrix operating in spin 
space. The Bogoliubov-de Gennes Hamiltonian ((4]) is the 
starting point when we in Sec. IIII B1 derive the appropri- 
ate reflection amplitudes for quasiparticles impinging on 
the superconductor interface. 



A. Matrix current 

We seek the right-going charge and spin currents in 
normal metal lead 1, and start by introducing the matrix 
current^ 



h. 



xfj (t) = 27rer^ (a\(t)a a {t) - SJ(i)6 a (*)) , (8) 



where e is the electronic charge, and t z is a Pauli matrix 
in electron-hole space: 



1 
-1 



(9) 
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Charge and spin currents are obtained from the matrix and with a vector of the Pauli matrices and their complex 
current ([5]) as follows: conjugates, as the diagonal elements. 

/cW=^(4a«W), (io) 



and 



I »(') = ^Z)Pafl<A,Pa(*)>, 



(11) 



respectively. Summations run over electron- hole, spin 
and mode space, and p is a matrix with diagonal struc- 
ture in electron- hole space: 



Pa/3 



r a/3 







For a slowly oscillating scatterer, it is conve- 
nient to express the scattering matrix in the Wigner 
representatio n 22 ! 33 ! 34 



t + t' 



(12) In this representation, the matrix current is: 
I 



(13) 



{hMt)) = 2^t. 



I JO „ J-oo Jo 



de± de^ 
2tt 



AW 



(14) 



The current is expressed in terms of the center and rel- 
ative time coordinates T = (t' + t")/2 and r = t" — t' , 
and the Fourier transform of the distribution function 

fl {r)= f^e-^f^e). (15) 

When the scattering matrix S(e; t) is a concatena- 
tion of multiple time-dependent scattering elements, the 
Wigner representation of S will also be an infinite sum of 
time and energy gradients^ 2 - The magnetization dynam- 
ics is slow as compared to the time an electron spends in 
the scattering region. In the adiabatic approximation, we 
assume the scattering matrix evolves on a much longer 
timescale than the typical dwell times of particles inside 
the scattering region. In this regime, we formally expand 
S asM 

S(e;t) = S (e:,t) + A{e;t) + O{d?S ) (16) 
I 



where So is the "frozen" or instantaneous scattering ma- 
trix, and the matrix A represents all first-order gradi- 
ent corrections to So resulting from the concatenation 
of time-dependent scattering elements that describe the 
device. Unitarity of S to all orders in time- and energy- 
gradients implies^ 

S A*+AS% = i (dtSod^ - d e S 8t$£) = \p [So; S\] , 

(17) 

where a Poisson bracket P{. \ .} has been defined to ease 
the notation. In the following, scattering matrix argu- 
ments (e; t) are omitted in places where there is no risk 
of confusion. 

To obtain a local (in time) expression for the matrix 
current (j 14[) . we Taylor expand S to first order in time 
derivatives, and obtain the matrix current 



(A,a/3W) - ^ T a/sX/ j (/ Q (e) - A( £ )) ^0,07^0,^7 + ^"7 ^0,^7 + Sq^A*^ - -P {So >ai ; S'o i/ 3 7 }^ 

+ \ (-0*A(e)) {So^dtSl^ - d t S , aj Slfr) \ + (d?S ) (18) 



where Eqs. (|16[) and JT7]) have been used. The matrix current in Eq. (1181) is exact to first order in frequency of 
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the pumping parameter. 

Finally, we observe that in the absence of a voltage 
bias, the gradient corrections to the frozen scattering ma- 
trix, represented by A, vanish from the matrix current. In 
electro-chemical equilibrium, when V = 0, f e (s) = fh{s), 
the first line of Eq. fp~8|) vanishes, and the pumped current 
is determined by the frozen scattering matrix. Naturally, 
the same is also true for the time-dependent theory based 
on Floquet scattering matrices. 27 



B. Scattering matrix for a N|F|S structure 

In this section, the scattering matrix formalism derived 
for N|S structures^ is applied to our N|F|S trilayer. As 
described in Sec. ITTJ the scattering description of a N|F|S 
structure is greatly simplified by inserting a fictitious nor- 
mal metal lead (N 2 ) between the two scattering regions, 
thereby spatially separating spin-dependent scattering in 
F and Andreev reflection at the N2IS interface. 29 The 
scattering matrix Sf, describing the disordered ferromag- 
netic region, is block-diagonal in electron-hole space. We 
write Sf as 



S F (e;t) 



s F (e;t) 







s F (-s;t)' 



where the diagonal elements are 



sf 



hi r 2 2 



(19) 



(20) 



Here, ru and iy are matrices in spin-space that describe 
reflection of an incoming electron in lead i, and transmis- 
sion of an electron from lead j to lead i, respectively. 

Electrons and holes with opposite spins are coupled 
by Andreev reflection at the superconductor interface, 
where an incoming electron (hole) is reflected as a hole 
(electron) with reversed spin direction. The reflection 
amplitudes are derived by matching propagating wave 
functions in N2 with evanescent wave functions in the su- 
perconductor. The resulting scattering matrix read a 29 ' 37 





he 







h _ 





-iaa y e~ 



\aa v e l << 




(21) 



where a = exp [— i arccos(e/A)]. 

The total scattering matrix of the N|F|S structure is a 
concatenation of Sf and r A , and in terms of the frozen 
scattering matrices, we obtain the familiar result a 29 ' 32 

S ee ( £ ;t) = r u ( £ ) 

+ t 12 (e)rf h (sy 22 (-e)M e (e)rt(e)h 1 (e), (22a) 
S$ h (e;t) = r* u (-e) 

+ tt 2 (-e)rt(e)r 2 2(e)M h (e)r A h (e)t* 21 (-e), (22b) 

S% h (s;t) = h 2 (e)M h (e)r? h (e)t* 21 (-s), (22c) 

S^(e;t) = tU-e)M e (s)ri(s)t 21 (e), (22d) 



where time arguments are omitted on the right hand side 
of the equations for sake of notation. Multiple reflections 
between S and F, mediated by propagations through N2, 
are described by 

M e {e) = [l-rt{e)r 22 {e)r* h {e)r* 22 {-e)]-\ (23) 
M h (e) = [l-r* h (ey 22 (-e)rt(e)r 22 (e)]-\ (24) 

From Eqs. (f2"2"j). and using rf h {—e)* = r^ e (s), one obtains 
the following symmetry relations for the total scattering 
matrix: 



and 



S ee (e;t)= [S hh (-e;t)}\ 



S £h (e;t)=[S he (-e;t)}' 



(25a) 



(25b) 



The frozen scattering matrices in Eqs. (]22j) are all time- 
dependent due to the slowly varying magnetization in 
the ferromagnet. Arguably the easiest way to evaluate 
the matrix current, is to perform a spinor rotation that 
aligns the spin quantization axis with the instantaneous 
magnetization direction! 16 i 21 The total scattering matrix 



Sb(e;t) 



qee qeh 

qhe qhh 
°0 °0 



(26) 



can be related to the total scattering matrix S_ in the 
rotating frame by the spinor rotations 



S (e;t) = W\t)S(e)W(t), 
where W{t) = V{t)U{t), with 



(27) 



_ (U(t) 



'exp [i^cr 2 ] 



exp [- 



and 



V[t) ~ { V(t) 



(28) 



exp 



2 u 



exp 



2 u 



(29) 

In the rotating frame, Sjj e and S_Q h are both diagonal in 
spin space, while S_Q h and Sjr] e , which mix spin a electrons 
with spin —a holes, only have off-diagonal elements. 

Now that the matrix current and relevant scattering 
matrices are derived, we proceed to study pumped charge 
and spin currents for a voltage biased trilayer structure. 



IV. PUMPED CURRENTS OUT OF 
EQUILIBRIUM 

A complication that arises when the system is driven 
out of equilibrium, is that time- and energy gradients 
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of the frozen scattering matrix must be evaluated. Be- 
fore presenting the detailed expressions for charge and 
spin currents in the normal metal lead, we derive the 
required gradient corrections. Due to electron- hole sym- 
metry (??), it is sufficient to consider only A he in the 
gradient correction. 



A. Gradient correction matrix 

In the following, we determine A he by a formal gradient 
expansion of the corresponding scattering matrix S he , 
whose full time and energy dependence of S he is given 
by (see Eq. (|22a|) ): 



where 



S he {e;t)= (tI 2 oM e or&ot 21 )( e; t). 



(30) 



Evaluating the convolutions in the Wigner represen- 
tation can be done by systematically expanding the 
exponentials^ 

(AoB)(e;t) = j( a ? B ?- B ? a ?)' a A{e;t)B( E ;t), (31) 

where the superscripts indicate which matrix the oper- 
ator works on. A significant simplification of the fi- 
nal result is achieved when the superconducting gap is 
much less than the exchange energy, A <C A xc ,Ep. The 
energy dependence is then only determined by the en- 
ergy dependence of the Andreev reflection. Since we are 
evaluating the energy gradients close to the Fermi level, 
d E sp -C d e r A , and we obtain the simplified expression 
for the gradient matrix A he : 



A = VUd t (m ■ er)W t V t = d t 9o* + sm9VLa y . 



(35) 



Multiplying r^ e , which is ~ a v , with A, and using that 
the other components in the equation are all diagonal, 
brings us to the conclusion that T e is diagonal in spin 
space. Finally, we note that T he — > for a vanishing 
ferromagnetic ordering parameter. 

Once the gradient corrections to the frozen scatter- 
ing matrix are derived, one can obtain non-equilibrium 
pumped currents to first order in pumping frequency. 



B. Pumped charge current 

According to Eq. ([ID]), the charge current is obtained 
by tracing the matrix current (|18p over electron-hole, 
spin and mode space. Making use of the electron-hole 
symmetries from Eqs. (I25a|) (I25b|) . and using that both 

Tr (d t S§ e So ef \ = and Tr {fit^S^} = 0, one finds 

that the pumped charge current is determined by 



Je([f e (e)-f h (s)]Tr{s^S^ 



(36) 



to first order in pumping parameter frequency. Using 
that A he = -id e d t Sft e + T he , the current (J3BJ) simplifies 
to 



A he {e-t) « --d £ d t S^ + it* 12 d £ (M e r£ e )d t t 21 

+ itl 2 d £ d t M e rtt 2 i + Yt* 12 d t M e d e M- l M e rttn 
-it* 12 M e rtd t r22d s r? h r* 22 M e rtt2i 

= - l -d e d t s* e + r he . (32) 

Here, Sq s is the frozen scattering matrix from Eq. (I22al) . 
and Before evaluating the currents, we observe that T he 
in the rotating frame is diagonal in spin space. This 
fact, which is important when evaluating non-equilibrium 
pumped charge and spin currents, can be seen from 



T he = UV f T he VLi, 



(33) 



with 



r e =2&de(M e ^ e )A(t 21t -t 2U ) 

+ ki 2 (Me t -M^)9e(^ e Ar 22 ri)r| 2 M e r4i 



;tl 2 M e rt(L 22t - L^kd^r^M^tt 



.21 



-^(Met-MeXeAfel, (34) 



IS) 



del [f e ( e )-f h (e)] 



e 

2^ 



x Tr{^ le ^ let + r he S% e1 + S£ e r he ^ } ) . (37) 



Any non-equilibrium pumped contributions to the cur- 
rent are determined by the remainder T he from Eq. (I34I) . 
However, as pointed out at the end of Sec. IIV Al T he is a 
diagonal matrix in spin space. From Eq. ()22a[) , we know 
that S_q & is strictly off-diagonal in spin space. This im- 
plies that Tr |r /le S , Q < ^j = 0, and the charge current is 
reduced to the stationary result: 



Ic = tt de[f e (e)-f h (e)]~g(e) 
where the total conductance is defined as 



(38) 



(39) 



The result in Eq. (|38|) shows that there is no pumped 
charge current in N|F|S structures, even when there is an 
additional bias voltage driving the system, e.g. there are 
no bilinear contributions proportional to the bias voltage 
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and the FMR frequency. The stationary result is similar 
to that obtained in FS|N structures^! a result that indi- 
cates that the total scattering matrix for a disordered re- 
gion coupled to a ferromagnetic superconductor, is struc- 
turally equivalent to that of a disordered ferromagnetic 
region coupled to a superconductor. The two structures 
have different scattering matrices, however, and therefore 
the expressions for the conductances differ. 



C. Pumped spin current 



We proceed by evaluating the pumped spin current to 
first order in pumping parameter frequency. Utilizing the 
electron-hole symmetry relations for the total scattering 
matrix, we obtain 



[ s(t) = ^J°°de (f e (e) - f h (e)) [Tr [a* (s' ie S^ + T he S^ + S^T he ^ } + 9 e ImTr {<r*d t S% e S^}' 

de (~dj e (e)) [imTr {ad t S% e S^} - ImTr {<T*d t S* e S^}] . (40) 



1 

47T 



Introducing the conductance polarization 



m n. v 



he 
—X.'X.mn 



ohe 

— ti,mn 



and the generalized mixing conductance^ 



~9 U 



nee nee* , qhe oh< 



(41) 



:*,»}. (42) 



we find the following expression for the spin current: 



Is(t) = 



i- f X As (Me) - f h (e))d e (m x d t m(g + 2Re £ SHl mn S^ mn ) + 29 t mlm^£^ mn S^ mi 

00 ^ m,n m.n 

1 f°° 

H / de d E f e (e) (m x (^mRec/^ + d t mlmg 

4tt 7.^ 



(43) 



The term ~ pgm(t) on the right hand side of Eq. (j43|) 
corresponds to the non-equilibrium bias voltage spin cur- 
rent observed also in the absence of a precessing magne- 
tization vector. Terms in the final line are similar to 
those derived previously within electro-chemical equilib- 
rium pumping theory for F|Ni&, and FS|N structures^. 
However, we ask the reader to note that the general- 
ized mixing conductance in Eq. (3) in Ref. [U is valid for 
triplet superconductors only; the correct mixing conduc- 
tance for a singlet superconductor is given by Eq. ((32) 
above. The remaining terms on the right hand side of 
Eq. (1431) are non-equilibrium, pumped contributions to 
the spin current. They depend on pumping parameter 
frequency via dtm and the A term from Eq. (|3"5j) . which 
is contained in the gradient remainder T he . 



V 



Finally, we would like to point out that there are no 
pumped contributions to the longitudinal spin current 
7a = m ■ I s . The terms in the second and third line of 
Eq. (|43| are transverse with respect to the magnetiza- 
tion m, so this leaves only a possible gradient remainder 
contribution coming from T he . However, due to the par- 
ticular matrix structure of T he mentioned in Sec. IIV Al 
m ■ Tr{<T*r ft-e S l Q e ^} vanishes. This observation implies 
that the longitudinal spin current is stationary and un- 
affected by the precessing magnetization. Thus, to first 
order in precession frequency: 

7]l = m(t)-I s (t) = --L J°° de(f e (e)-f h (e))p~g. (44) 

In the following, we will investigate pumped charge and 
spin currents when the ferromagnetic region is longer 



8 



than the typical transverse spin coherence length. 



D. Long ferromagnet limit 

When the length Lp of the ferromagnet is longer than 
the transverse spin coherence length, 



•10" 



Lp > Xi 



k 



Fl 



(45) 



where kp a is the Fermi wave vector of a spin a electron, 
we expect to find a mixing conductance that is deter- 
mined by the properties of the N-F subsystem, charac- 
terized by the spin-dependent conductances^ 



m.n 



(46) 



Indeed, in the limit (|45|). one can disregard "mixing 
transmission" terms, J2 m ,n i<r,mrd£-<r,mn 0' so tnat 
T, m ,nS!llmn^l* mn 0. Disregarding interference 
terms between reflected and transmitted electronic wave 
functions, one obtains 



Ertee qee* \ " 



If, mn— mm 



(47) 



for a long ferromagnet. This implies that — > g^, 
while the total conductance g and the conductance po- 
larization p remain unchanged. Since the mixing con- 
ductance is now determined by properties of the N-F 
structure, energy gradients of the mixing conductance 
should be disregarded in the limit A <C A xc ,Ep, as de- 
scribed in Sec. IIV Al Finally, by an explicit calculation, 
one can show that Tr{T he a*} ~ t a t*_ a , which van- 
ishes when Eq. (|45[) holds. To summarize, when the fer- 
romagnet is longer than the transverse spin coherence 
length, the charge current and longitudinal spin current 
are still given by 



e 

2^ 



de(/ e (e) - fh(e))g, 



(48) 



and 



jii = -JL 

An 



de(f e (e) - f h (e))pg, (49) 

while the transverse spin current is simplified to 
1 f°° 

I s( t ) = -g- / de(d e f e (e) ~ dj h (e))gm x d t m 
1 f°° 

+ — / ded e f e (e) (Reg ti m x d t m + Jmg'^d t m) . 

J — oo 

(50) 

With no applied bias voltage, the pumped spin current 
in Eq. (|50p is identical to that found in N-F systems^, as 
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(a) Transverse spin component dtm ■ I s /Cl 
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Lf/^f 

(b) Transverse spin component (m X dtm) ■ I s /Cl 

FIG. 2: Exact (red line) and approximate (blue dashed line) 
transverse spin currents for a ballistic N-F-S structure as 
functions of length of the ferromagnetic region. In the plot, 
E F = 10 eV, A xc = 9E F /16, A = E F /160, eV = A/2 and 
Q = 0.2 GHz. 



should be expected. In this situation, emission of spins 
from the ferromagnet into the normal metal are unaf- 
fected by the superconductor. 

To compare the exact result (|4"5)) with the long ferro- 
magnet approximation of Eq. (|5T)|) . we plot in Fig. [5] the 
spin current along dtm for a ballistic N|F|S trilayer, as 
a function of the ratio between the ferromagnet length 
(Lp) and the transverse spin coherence length (Xf) de- 
fined in Eq. (1431) . When Lp < Xp, non- negligible "mixing 
transmission" terms combine with energy gradients of the 
scattering matrix and produce large deviations between 
the two equations. As exceeds Xf, the fit improves 
and the exact result oscillates towards the spin current 
obtained by the approximate Eq. (|50[) . 
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V. CONCLUSION 

In conclusion, we have derived non-equilibrium 
pumped charge and spin currents to first order in pump 
frequency, using time-dependent scattering theory. Mag- 
netization precession induces transverse spin currents, 
but neither charge nor longitudinal spin currents, which 
are both given by their stationary values. The currents 
are expressed in terms of generalized, spin dependent con- 
ductances, that include spin-dependent scattering in the 



ferromagnet and Andreev reflection at the F|S interface. 
Finally, we consider trilayers where the ferromagnetic re- 
gion is longer than the transverse spin coherence length, 
and derive an approximate expression for the transverse 
spin current. Numerical calculation of the spin current in 
a ballistic trilayer shows good agreement between exact 
and approximate spin currents for ferromagnets whose 
layer thicknesses exceed the transverse spin coherence 
length. 
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